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> . 1 INTRODUCTION 
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The integrated Sachs-Wolfe (iSW) e ffect jSachs & Wol"fej|l967l ; 
iHu & Sugivamal 1 19941 ; ICooravl I2002T) . which refers to the fre- 
quency change of cosmic microwave background (CMB) pho- 
tons if they cross time evolving gravitational potentials, is a 
direct probe of dark energy because it vani shes in cosmolo- 
gies with £l m = 1 dCrittenden & Turokl 1 19961) . By now, it has 
been detected wi th high significance with a number of differ- 
ent tracer objects iFosalba et alj2003l:lBoughn & Crittendenll 20041: 



m 

Nolta et alj 20041: IPa dmanabha n et al 
20061: ' 



2005 



Giannantonio et al. 

Pietrobon et al. 2006; Gaztanaga et al. 2006; Cabre et al. 
Vielva et a l. 2006; Rassat et al. 20071 iMcEwen et alj|2007l ; 
Giannantonio et aul2008h , and derived parameter constraints pro- 



2006; 



ABSTRACT 

The subject of this paper is an investigation of the nonlinear contributions to the spectrum 
of the integrated Sachs-Wolfe (iSW) effect. We derive the corrections to the iSW-auto spec- 
trum and the iSW-tracer cross-spectrum consistently to third order in perturbation theory and 
analyse the cumulative signal-to-noise ratio for a cross-correlation between the PLANCK and 
EUCLID data sets as a function of multipole order. We quantify the parameter sensitivity and 
the statistical error bounds on the cosmological parameters Q m , <x 8 , h, n s and w from the lin- 
ear iSW-effect and the systematical parameter estimation bias due to the nonlinear corrections 
in a Fisher-formalism, analysing the error budget in its dependence on multipole order. Our 
results include: (z) the spectrum of the nonlinear iSW-effect can be measured with 0.8cr sta- 
tistical significance, (if) nonlinear corrections dominate the spectrum starting from t 10 2 , 
(Hi) an anticorrelation of the CMB temperature with tracer density on high multipoles in the 
nonlinear regime, (;V) a much weaker dependence of the nonlinear effect on the dark energy 
model compared to the linear iSW-effect, (v) parameter estimation biases amount to less than 
O.lcr and weaker than other systematics. 

Key words: cosmology: CMB, large-scale structure, methods: analytical 



ear iSW-effect, the RS-effect shows a flatter spectral depen- 
dence and dominates the signal at higher multipoles exceeding 
I i 100. Analytical, perturbative derivations agree well with the re- 
sults from n-body simulations jTuluie et alj[l996l : ICai et alj|2008l : 
ISmith et al.1 120091 : ICai et alj l2010h . The~non-Gaussianities intro- 
duced into the CMB by the nonlinear RS-effect are very weak 
jMollerach elat] Il995l: llvlunshi et al.1 Il995l ; ISpergel & Goldberd 
ll999l : lGoldberg & Spergell 19991 although the first two papers work 
in the context of a SCDM-cosmology, their results are still ap- 
plicable to ACDM). The RS-effect from the local Universe has 
been found to am ount to ~ 2ii.K in the most massive structures 
dMaturi et alj|2007h forming in a constraint realisation. 



vide support for a ACDM cosmology. 

C ontrarily, the nonlinear iSW-effect , or Rees-Sciama (RS) 
effect teees& Sciamalll968l : ISeliaklll996l : ISchafer & Bartelmannl 
2006) is difficult to detect and shows only a w eak signal amout- 
ing to < 2cr i n the spectrum JCooravl I2002T) or up to 0.8cr in 
the bispectrum jSchaferll200i) . The cross-correlation with weak 
lensing has been shown to be fe asible, but weak with cur- 
rent surveys jNishizawa et alj l2008h - In comparison to the lin- 
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The topic of this paper is the contamination of the iSW- 
spectrum by the nonlinear RS-spectrum at intermediate multipoles: 
In a measurement of the linear iSW-effect, nonlinear contributions 
will alter the shape of the observed spectrum and can affect the 
estimation of cosmological parameters by introducing estimation 
biases. We investigate dependence of parameter accuracy as well 
as the parameter estimation bias as a function of maximum mul- 
tipole order considered. Specifically, we use a Fisher-matrix ap- 
proach to quantify the statistical and systematical errors, analyse 
the error budget as a function of multipole order and derive the 
optimal maximum multipole moment which minimises the com- 
bined error for individual parameters. The nonlinear iSW-effect is 
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the most important contaminant at intermediate multipoles, with 
the kinetic Sunyaev-Zel'dovich effect starting to dominate at higher 
multipoles above thousand. 

After summarising key formulae describing structure forma- 
tion in dark energy cosmologies in Sect. [2] we introduce line of 
sight expressions of the two relevant observables in Sect. [3] We 
carry out a perturbative expansion of the source fields to third order 
in Sect.|4]and derive the spectrum C Ty (t) between iSW-temperature 
perturbation r and the galaxy density y to third order in Sect.[5] We 
quantify the degeneracies between the cosmological parameters us- 
ing a Fisher-matrix analysis in Sect.|6]and extend this formalism to 
describe the parameter estimation bias in Sect. [7] A summary of 
our results is compiled in Sect. [8] 

As cosmologies, we consider spatially flat homogeneous dark 
energy models with constant dark energy equation of state, and 
with Gaussian adiabatic initial conditions in the cold dark matter 
field. Specific parameter choices for the wCDM-fiducial model in 
the Fisher-matrix analysis are H = 100/7km/.?/Mpc with h = 0.72, 
n„, = 0.25, a b = 0.04, erg = 0.8, w = -0.9 and n s = 1, with con- 
stant unit bias for the tracer galaxy population. 



2 COSMOLOGY AND STRUCTURE FORMATION 
2.1 Dark energy cosmologies 

In a spatially flat dark energy cosmology with a constant dark en- 
ergy equation of state parameter w, the Hubble function H(a) = 
dlna/d/ is given by 



H 2 (a) £1„, 1 - Sl m 
— + 



7 3(l+u>) 



(1) 



The value w = -1 corresponds to the cosmological constant A. 
The conformal time, which is related to the cosmic time t by the 
differential d^ = At I a, follows directly from the definition of the 
Hubble function, 



n = f da — , 



(2) 



in units of the Hubble time t H = 1 /H . Correspondingly, the defini- 
tion of the comoving distance x ls given by x — CJ 1 with the speed 
of light c. 

2.2 CDM power spectrum 

A common parameterisation for the CDM power spectrum is 
P(k) oc k"'T 2 (k) for describing the Gaussian fluctuation statistics 
of the homogeneous and isotropic cosmic density field 6, 

(S(k)S(k'y) = {2nfS D {k - k')P(k) (3) 

According to Bard een et al] (l986), a convenient fit to the CDM 
transfer function T(k) is 

ln(l + 23Aq) 



T{q) 



2.34a 



■ (l + 3.89a + (16.1a) 2 + (5.46a) 3 + (6.71a) 4 ) 5 



where the wave vector q is given in units of the shape parameter T =; 
£l,„h. P(k) is normalised to the value cr 8 on the scale R = 8 Mpc/A, 



t2r = ihl Ak ^ w ' {kR "> p ^ 



(4) 



with a Fourier-transformed spherical top-hat W(x) = 3ji(x)/x as 
the filter function. j ( (x) denotes the spherical Besse l function of 
the first kind of order t jAbramowitz & Stegur]|l972h . [Srnith et al.l 



j2009h found that nonlinear effects in the biasing model amount to 
~ 10%, but for simplicity, we assume a linear, local, non-evolving 
and scale-independent biasing scheme, 



A/7 



Ap 
P ' 



(5) 



and relate fluctuations An in the spatial number density n of galax- 
ies directly to the dark matter overdensity S = Ap/p. 

2.3 Structure growth in dark energy cosmologies 

The linearised structure formation equations, i.e. the continu- 
ity, Jeans and Poisson equations, can be combined to the 
growth equation (jWang & Steinhardtll 19981 ; iTurner & Whitell 19971 : 
iLinder & Jenkinsl2003l) ,~ 



d 2 1 / dln#\ d 

-rD + + - 3 + — D + 

da 1 a\ din a / da 



2a- 



Q m (a)D+(a). 



(6) 



whose solution D + (a) describes the homogeneous growth of the 
density field, S(x, a) = D + (a)6(x, 1). In the standard cold dark mat- 
ter (SCDM) cosmology with fi„, = 1 and 3 + d In H/d In a = | , this 
solution is easily derived to be D + (a) = a. This motivates the choice 
D + (0) = and d/daD + (0) = 1 for the initial conditions, due to 
matter domination at early times. The second solution D_(ci) = 1/a 
decays rapidly and has no influence on the late-time iSW-effect. 



3 OBSERVABLES: ISW-EFFECT AND TRACERS 
3.1 iSW-temperature perturbation 

The iSW-effect is caused by gravitational interactions of CMB pho- 
tons with time-evolving potentials O. The fr actional perturbation r 
of the CMB tempera ture 7cmb is given by dSachs & Wolfein"967l : 
iRees & Sciam3l 19681) 



AT 



2 C XH , ao 
-- dxa 2 H(a) — , 
c 3 Jo da 



(7) 



The gravitational potential <I> is a solution to the comoving Poisson 
equation, 



AO = — °—S, 
2a 



(8) 



Substituting into the line of sight expression for the linear iSW- 
effect t (integrating along a straight line and using the flat-sky ap- 
proximation) yields 



3Q,„ 



C XH 9 d /D + \ 



Xh 



(9) 



c Jo 

where the inverse Laplace operator A~' lx 2 H solves for the potential: 

A- 1 

v = —6. (10) 

Xh 

The square of the Hubble distance xh = c/Ho makes the differential 
operator dimensionless. 



3.2 Galaxy density as a large-scale structure tracer 

The projected galaxy density y can be related to the CDM density 
5 via 



d XP (z) — D + (x)S, 
Jo &X 



(11) 
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where p(z)dz is the redshift distribution of the surveyed galaxy 
sample, rewritten in terms of the comoving distance x- We use 
the redshift distribution of the main galax y sample of EUCLID 
dRefregier & the DUNE collaboration 2008), which will / skv = 0.5 
of th e sky with a median redshift of Zmed = 0.9 JDouspis et al.l 
2008). We use the parameterisation proposed bv lSmail et alj fl995) 



p(z)dz = po 



-,exp 



— ] I dz with — 

zo 



1 

Po 



z 4r 3 - 



(12) 



for p(z)dz, with zo = 0.64. We assume a constant bias of b = 1, 
which we absorb into the normalisation cr 8 of the power spectra. 



4 PERTURBATIVE CORRECTIONS 

For a consistent derivation of the iSW-spectrum including correc- 
tions due to the nonlinearly evolved source fields one needs to carry 
out a perturbative expansion to third order, 



6{x, a) =s ^ DJCfiO^Cx) + 0(6 4 ). 



(13) 



The linearity of the Poisson equation conserves the perturbative se- 
ries, 



' da a 



(14) 



and suggests that the time derivative of the potential is oc 
d(D" /a) /da. In perturbation theory, the second and third order cor- 
rections to the density field are given by 



V_ 



M 2 (k-p,p)S(p)S(\k-p\), 



(15) 



J (2n 

<* <3) (*) = f 7S7 f ^M 3 (p,q,k-p-q)d(p)S(q)6(\k-p-q\), (16) 

where the mode coupling functions M 2 (p, q ) and M 3 (p,q,r) (see 
ISahni & Colesll 19951 : Bernardea u et all 20021) are a consequence of 
the inhomogeneous growth and introduce non-Gaussianities in the 



evolved density field. The power spectrum PgP(Ji) 
density field thus acquires the corrections 



P(k) of the 



1 f P^M 1 (k-p,pfP(\k-p\)P(\p\), 

1 f ll-M 3 (k,p,-p)P(k)PQp[y 
J (2tt) 3 



(17) 



(18) 



In the computation of these corrections, the cylindrical symmetry 
of the kernels Mo and M3 can be taken advantage of, reducing to 
a twofold integration 2np 2 dpdp with p being the cosine of the an- 
gle between k and p. Fig.[TJshows the time evolution of the source 
fields, i.e. the growth function D" + (a) for the density field, and the 
time derivative of D" + (a)/a for the iSW-source field, both up to per- 
turbative order n = 3. While the growth functions show a similar 
behaviour in higher order, the derivatives are qualitatively very dif- 
ferent. The evalutations of the integrals is done in a coordinate sys- 
tem whose p,-axis is parallel to k z . The nonlinear corrections to the 
CDM spectrum P(k) are shown in Fig. lAll 

An interesting peculiarity of the nonlinear RS-effect in com- 
parison to the linear iSW-effect is worth mentioning: Whereas in 
SCDM-cosmologies the iSW-effect vanishes due to D + (a) = a and 
is nonzero in dark energy cosmologies, the RS-effect is strongest 
in SCDM and weaker in dark energy cosmologies, at least at the 
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Figure 1. Time evolution of the source term D" + (a) for the density field 
(thick line) and the modulus of d(D" /a) /da for the iSW-effect (thin line), 
for the linear order n = 1 (solid line) and the nonlinear corrections n = 2 
(dashed line) and n = 3 (dash-dotted line), with ACDM as the cosmological 
model. 



low redshifts we observe, where D + (a) =2 a" with a < 1. Further- 
more, the cross-spectra of the RS-effect are proportional to £l m a\ 
(up to third order in perturbation theory), in contrast to the iSW- 
spectrum, which scales as Q. m a\. The dependence on the dark en- 
ergy eos-parameter w is weaker in the nonlinear effect and the 
shape of the spectrum (determined by h and n s ) becomes less im- 
portant because of the integrations over d 3 p carried out in pertur- 
bation theory. These arguments motivate the quantification of the 
RS-contamination of the iSW-spectrum, and their interference with 
the estimation of cosmological parameters. 



5 ANGULAR POWER SPECTRA 

In summary, the line of sight integrals for the iSW-temperature per- 
turbation r (n) and the galaxy density y <n) in order n read: 



d X a 2 H(a) 



Jo 

d XP {z)-^Dl(x)& n \ 
d X 



c 

rXH 



(19) 
(20) 



where we have defined the dimensionless potential <// n) = 
A~ 1 6 m /Xh fr° m the inversion of the Poisson equation, rescaled 
with the square of the Hubble distance xh = c/Hq for convenience. 
Due to the linearity of the Newtonian Poisson-relation, the pertur- 
bative corrections in 6 map directly onto the corrections in ip. The 
weighting functions 



3n 



d D", 



da a 



dx 



(21) 
(22) 



can be identified, which allow the expressions for the angular cross 
spectra to be written in a compact notation, applying a Limber- 
projection jLimbeJI 19541) in the flat-sky approximation. 



00 



r 

Jo 



% W^\ X )W^(x)P^\l/x) 
X 



(23) 
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aw® = r % WOW 

Jo X 



*J H*l 

4 w< 2) Cr)M4 2) 0r)^ 2) (^) (25) 
o X 

with the cross-spectrum P^?(k) = Pgg(k)/(xHk) 2 . The expression 
for the spectrum C^XO has been symmetrised. The angular auto- 
spectra of the temperature perturbation r are given by: 

(26) 
r*n Ay 

C£ 3) (€) = 2 J o -£ ^(tfW^fr^W*) ( 27 ) 
CgO W = jf** 4| W?(x?^\ei X ) (28) 

In analogy to F$(K), the spectrum of the potential <p is defined as 

P^ik) = P%\k)l(XHkf. Finally, the spectra of the galaxy density 
y can be evaluated to be: 

h, X 

rXH 



Jo 



J"XH Ay 
o A" 



r%w?(x?p%xm 

Jo X 



(29) 
(30) 
(31) 



Collecting all terms, the full spectra consists of one first and two 
second order contributions, 



C Ty {£) 

C TT {1) 

c„(t) 



C[\ v \£) + C< 22) (f) + C { T f{€), 

c^\e) + c< 2 T 2) (fi + c^Xe), 
c^Xe) + c^XO + c^Xo. 



(32) 
(33) 
(34) 



Figs. [2] and [3] give the iSW-auto and cross-spectra, respec- 
tively, split up into linear contributions and the two perturbative 
corrections. The iSW-auto spectrum is dominated on multipoles 
larger than 100 and the cross-spectrum is suppressed by the neg- 
ative correlation between iSW-effect and tracer density on sim- 
ilar scales, leading to a sign change of the cross-spectrum at 
£ - 500, which confirms earlier perturbative and w-body results 
(Selia l3ll996l:ICooravll2002l : lNishizawa et alj2008tlCai et alj|2008l : 
ISmithetal.ll2009l) . but using a different perturbation theory ap- 
proach. arsinh(jr) is equal to x for \x\ <k 1 and oc In x for \x\ » 1, 
which allows to show the logarithmic behaviour of C Ty {£) despite 
the sign change. Another interesting feature is the fact that the non- 
linear effect is much less sensitive on the choice of cosmological 
parameters, in particular the dark energy equation of state param- 
eter w, for which the RS-spectra depicted differ by about 5%. The 
sign change and its sensitivity on w is mostly driven by changes 
in C^y (€). Fig. IA2I gives the cross-spectrum in a logarithmic rep- 
resentation and shows that the anti-correlation between t and y is 
a generic feature of nonlinearly evolving structures from angular 
scales of i ^ 70 on, but the linear effect shifts the anticorrelation 
scales to much higher multipole moments. The sign change can 
be easily explained by the fact that in linear structure formation 
potentials are constant or decay slowly, depending on cosmology, 
whereas in nonlinear structure formation the potentials grow fast, 
which manifests itself in the iSW-effect by causing temperature per- 
turbations of opposite sign. 




multipole order I 

Figure 2. Angular iSW-spectrum C TT (() of the iSW-effect (solid line), 
split up into the linear effect cS4''(0 (dashed line) and the nonlinear RS- 



corrections C^Xt) + CYr'if ) (dash-dotted line). The plot compares spectra 



-(13) 



for vi'CDM with w 
lines). 



-0.9 (thick lines) with ACDM with w 



-1 (thin 
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Figure 3. Angular iSW-cross spectrum C Ty (£) of the iSW-effect (solid line), 
split up into the linear effect Qi'(f) (dashed line) and the nonlinear RS- 
corrections C$. 22) (f)+C^, 3) (£) (dash-dotted line), both for wCDM (w = -0.9, 
thick lines) and ACDM (w = -1, thin lines). 



6 STATISTICAL ERRORS 

In this chapter, we recapitulate the estimation of statistical preci- 
sion on p arameters derived fro m angular iSW-spectra with Fisher 
matrices iTegma rk et alj[l997h . and the accurac y of the parameter 
estimation with an extended Fisher formalis m dCabre et alj|2007l : 
lAmara & RefregieifeOollTaburet et alJl2009h . 



6.1 Fisher-matrix for the iSW-spectrum C Ty (^) 

The Fisher matrix, which quantifies the decrease in likelihood if 
a model parameter x tl moves away from the fiducial value, can be 
computed for a local Gaussian approximation to likelihood £, oc 
exp(-^ 2 /2). The Fisher-matrix for the measurement of C Ty (€) is 
given by 
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dx„ 



(c ry (/),c, 7 (o) 



dC Ty (£) 

dXy 



(35) 



We construct the Fisher-matrix F^, for ACDM as the fiducial cos- 
mological model, with fiducial values for the parameters being 
n„, = 0.25, <Tg = 0.8, h = 0.72, n s = 1 and w = — 1. Implicitly, 
we assume priors on spatial flatness, £2„, + fl A = 1 and neglect the 
weak dependence of the shape parameter on the baryon density Q. b . 
CMB-priors on the cosmological parameters are incorporated by 
adding the CMB Fisher matrix F™ B , 



(36) 



6.2 Noise modelling 

In an actual observation, the iSW-power spectrum is modified by 
the intrinsic CMB-fluctuations, the instrumental noise and the beam 
as noise sources, assuming mutual uncorrelatedness of the individ- 
ual contributions. The galaxy correlation function assumes a Pois- 
sonian noise term, 



'ttCO 
yy(t) 



CrAO + CcmbW + W- T [ B- 
C yy {() + -, 

n 



(37) 
(38) 



For PLANCK's noise levels the value w? 1 = (0.02pK) 2 has 
been used, and the beam was assumed to be Gaussian, B~ 2 (t) = 
exp(A6» 2 t(l + 1)), with a FWHM-width of AO = I'.l, correspond- 
ing to channels of PLANCK closest to the CMB-maximum at 
=: 160 GHz. 

EUCLID is designed to survey the entire extragalactic sky and 



to cover the solid angle Afl = 2tt, corresponding to f s 



sky 



0.5, 



yielding a total of n = 4.7 x 10 galaxies per steradian at a density 
of 40 galaxies per squared arcminute. The observed cross power 
spectra are unbiased estimates of the actual spectra, 



CUt) = CUt), 



(39) 



in the case of uncorrected noise terms. We determine the spectrum 
Ccmb(I) °f me pr imary CMB anisotropies with the CAMB code 
dLewisetal., 2000). The covariance of the spectrum C Ty (t) is given 
in terms of the observed spectra C TT (£), C Ty (£) and C yy (l) which 
follow directly from applying the Wick-theorem, 



Cov(C ry ,C ry ) = j- [Cly{€) + C TT ({)C„({)] 



(40) 



In all applications considered in this paper, PLANCK causes the 
dominating noise contribution in comparison to the Poisson noise 
in the galaxy number density given by EUCLID. 

6.3 Detectability of the RS-effect 

The signal to noise ratio S of the cross-spectrum C Ty {t) reads: 
CUt) 



Cov(C rr (f))' 



(41) 



for mutually uncorrelated modes as in the case of a full-sky obser- 
vation. Fig. [4] shows the signal to noise ratio £ of a measurement 
of C ry (f ) including the nonlinear contribution at high I. The fig- 
ure suggests thay ideal cosmic variance limited experiments can in 
fact detect the RS-effect with a significance of 3.22<r (correspond- 
ing to a confidence of 0.998) integrating over all multipoles up to 
I = 3x 10 3 , and that this significance is reduced by the finite resolu- 
tion and the noise of PLANCK to a mere 0.77cr (0.558 confidence). 
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Figure 4. Signal to noise ratio of a measurement of the cross-spectrum 
C Ty (t), cumulative Y. (thick lines) and differential dS/d^ (thin lines), for 
the linear iSW-effect (dashed line) and the nonlinear RS-eftect (solid line). 
The plot compares the signal to noise ratio attainable with hypothetical 
cosmic-variance limited experiments (upper set of lines) with that reach- 
able by combining PLANCK with EUCLID (lower set of lines). 



Thus, the signal to noise ratio of the nonlinear effect is roughly 
smaller by an order of magnitude compared to that of the linear 
iSW-effect. Apart from the increasing correlation noise at high t it 
is the smallness of the spectrum around the cross-over scale which 
does not provide enough signal for a detection. Between t = 30 and 
I = 100 the cumulative signal to noise ratio sta gnates, which is not 
included in the computation bv lCooravl 120021) as his perturbative 
approach is not able to reproduce the small values of dE/d/ due to 
the sign-change of C Ty {€) at £ ^ 70. Despite the sensitivity of the 
cross-over scale on e.g. the dark energy equation of state parameter 
w it would be very difficult to measure this scale as the signal to 
noise ratio of each multipole is about 10~ 3 and as the same knowl- 
edge on w can be already derived from much smaller multipoles 
with sufficient accuracy. 



6.4 Parameter bounds and degeneracies 

The ^-function for a pair of parameters (x^,x v ) can be computed 
from the inverse (F~ v )u V of the Fisher matrix, 



Ax v 



(F-% (F~\ 



A.Xy 



(42) 



where Ax„ =x M - rf CDM . 



y/(F-%,(F- l ) v - 



The correlation coefficient r, n . is defined 



(43) 



and describes the degree of dependence between the parameters 
and x v by assuming numerical values close to for indepen- 
dent, and close to unity for strongly dependent parameters. The 
degeneracies between the cosmological parameters Q,„, cr 8 , h, n s 
and w estimated from the linear iSW-effect is shown in Fig. [5] 
along with the correlation coefficient from a measurement com- 
bining PLANCK and EUCLID data up to very high multipoles of 
i = 3000, including a prior from CMB data. 
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Figure 5. Constraints on the parameters fi,„, erg, n s , h and w, from the cross correlation of PLANCK with EUCLID. The ellipses correspond to 1 . . . 5<r 
confidence regions. Additionally, the vectors (S^,S V ) indicate the bias in the estimation of the cosmological parameters due to the nonlinear contributions and 
have been enlarged by a factor of 10. The estimation bias was derived for a multipole range extending up to £ max = 3 X 10 3 . The number in the upper right 
corner of each panel gives the correlation coefficient r™. The constraints include a prior from the statistics CcmbC) of primary CMB temperature anisotropies 
on both the statistical and systematical error. 



7 SYSTEMATICAL ERRORS 

In this section, we quantify how the interpretation of the data with 
the pure iSW-spectrum affects the estimation of cosmological pa- 
rameters, if in reality there are nonlinear RS-contributions at higher 
multipoles. Using this formalism, we seek to minimise the com- 
bined statistical error by finding an optimal angular scale £ opt down 
to which the iSW-measurement should be carried out. The nonlin- 
ear iSW-effect is the dominant contamination of the iSW-spectrum 
at intermediate multipoles, with the kinetic Sunyaev-Zel'dovich ef- 
fect becoming important at multipoles above I = 10 3 . The parame- 
ter estimation bias formalism has been validated with Monte-Carlo 
Markov-chains an d was found to be an excellent approximation for 
weak systematics faburet et all20ld) . 



7.1 Estimation bias formalism 

The angular iSW-spectrum C Ty {£) = Cisw(^) + CrsOT can be sepa- 
rated into the linear part dswiO and an additive systematic Crs(^) 
due to the nonlinear corrections, 



(44) 



CrsW = C<f (0 + c£ 3 >(0 (45) 

Using these relations, we define the power spectrum of the true 
model C,(() including nonlinear corrections, 

Cm = C.swW + CrsW = C^Xt) + C<f (f) + C<f (/) (46) 

as well as the spectrum of the false model Cf(€), which neglects 
these RS-contributions, 



c f (£) = CbwW = c^Xe), 



(47) 



where the observed spectra €/(£) are unbiased estimators of the 
theoretical spectra Cj(£) in each case, because of uncorrelated er- 
rors in each observational channel in the cross-correlation measure- 
ment method. The estimation of cosmological parameters is carried 
out from maximisation of the ^ 2 -functionals of the two competing 
models, 



w (cm -cm) 
xt = Zj 



cov [cm,cmv 
w [cm-Cfd)) 2 
4i covfc/w.c/co] 



and 



(48) 



(49) 
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i.e. the data is in reality described by C t (l) and, in the second case, 
fitted wrongly with Cf(£) instead of C,(£). The best-fit parameters x 
for each model can be derived by solving the equations (d^/dx^) = 
following from the respective ^-functional. 

For deriving the distance Xf — x t between the best-fit values 
of the true and the false model, we exp and the x\ function at the 
best-fit position x, in a Taylor series (see lTaburet et all2009h 



x )(x f ) = x )(x,) + £ JLftx,) 6, + \ £ q^X 2 M<) S M 6 r , (50) 



where the parameter estimation bias vector S =Xf—x, was defined. 
The best-fit position Xj of x\ can be recovered by extremisation of 
the ensemble-averaged (Xj)> yielding 

a 2 



dx^dx. 



which is a linear system of equations of the form 



Z G " 



(51) 



(52) 



where the two quantities G^, and follow from the derivatives of 
the ^-function, evaluated at x t , 



2 Gov- 



<max 

I 



- LRs(t) 



Cov~ 



dx^ dx v 



dx^dxy 



C RS (0 



dx u 



(53) 



The CMB priors can be incorporated by adding the Fisher-matrix 
F^T to G,,, 



Guv ■ 



<V +1%, , 



(54) 



for independent iSW- and CMB-likelihoods. The biases in parame- 
ter estimation from the iSW-effect are depicted alongside the de- 
generacies in Fig. [5] for a maximum multipole order of £ max = 
3000. The parameter estimation biases in the combined set of cos- 
mological parameters are very small, due to the weakness of the 
RS-signal in comparison to that of the iSW-effect, and due to 
the strong prior from primary CMB fluctuations. Typical values 
for misestimates in cosmological parameters are of the order of 
< O.lcr, and are negligible in comparison to statistical errors. 



7.2 Contamination of the iSW-spectrum 

In this section we consider the application of the iSW-effect for pro- 
viding independent constraints on individual cosmological param- 
eters. If the iSW-likelihood is combined with the CMB-likelihood 
according to eqn. {54}, the latter is by far dominating due to larger 
signal to noise ratio. Although the signal strength of the iSW-effect 
is not enough for fully constraining a standard dark energy cos- 
mology with five or more parameters, it is sufficient to place com- 
petitive bounds on single cosmological parameters. Therefore, we 



define the conditional systematical error cr^ 



and the sys- 



tematical error on a single parameter x^ while all other parame- 
ters are assumed to coincide exactly with their fiducial values. 

At low multipoles f max the error budget will be dominated 
by statistics, while the systematics due to the nonlinear contribu- 
tions are negligible. Conversely, the extention of the computation 
to higher multipoles i max will reduce the statistical error, but the 
RS-contributions will start to deteriorate the parameter accuracy. 




10 30 100 300 

multipole order ( 



1000 3000 



Figure 6. Conditional statistical errors <7),(dax) (thick lines) and systemat- 
ical errors bMmsa) (thin lines) for individual cosmological parameters: £l„, 
(solid line), <xg (dashed line), h (dash-dotted line) and w (dotted line) as a 
function of maximum multipole order ^ max . No prior information is used 
in the computation of the errors. The plot compares a cosmic variance lim- 
ited measurement (upper lines) with the cross correlation of PLANCK with 
EUCLID (lower lines). 



Fig.[6]depicts the individual conditional statistical and systematical 
errors as a function of maximum multipole order £ max . In compar- 
ison to the statistical errors on parameters derived with the iSW- 
spectrum, which are monotonically decreasing, the parameter es- 
timation biases due to RS-contributions have a more complicated 
behaviour with multipole order i, but remain always small in com- 
parison to the statistical error by more than one order of magnitude, 
for both cosmic variance dominated experiments and the combina- 
tion of PLANCK with EUCLID. The worst case is the constraint 
on w in a cosmic variance limited experiment, where the systematic 
error amounts to 20% of that of the statistical error. There are cer- 
tain scales at which the systematical errors are very small, namely 
as they change their signs, in agreement with changing parameter 
degeneracies on different angular scales. 



8 SUMMARY 

The topic of this paper is an investigation of the contamination due 
to nonlinearly evolving structures on the linear iSW-effect, and the 
consequent parameter estimation biases. 

(i) The angular spectrum of the Rees-Sciama effect was com- 
puted in third order perturbation theory. The spectrum C TT (€) of 
the RS-effect starts dominating that of the linear integrated Sachs- 
Wolfe effect from multipoles of t - 100 on. In particular the cross- 
spectrum C Ty (t) shows a sign change suggesting that the CMB 
temperature is anticorrelated with the galaxy density on nonlinear 
scales. This scale bears some sensitivity on w and is shifted from 
i - 70, where the sign change occurs in the spectrum of the nonlin- 
ear effect to i =i 500 for the combination of the linear and nonlinear 
effect. The sensitivity of the RS-effect on the dark energy model it- 
self is rather weak. 

(ii) By combining the PLANCK and EUCLID data sets one can 
measure the nonlinear RS-effect with a significance of 0.77cr out 
to multipoles of I = 3000, where the most important limitations 
are cosmic variance and PLANCK's instrumental noise. An ideal 
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experiment only subjected to cosmic variance would be able to de- 
tect the effect with 3.22cr, but at higher multipoles, confusion with 
the kinetic Sunyaev-Zel'dovich effect would occur. Measurements 
of the angular scale on which the sign change occurs are almost 
impossible as the signal to noise ratio on these scales is m 10~ 3 . In 
summary, the significance of the RS-effect is smaller to that of the 
iSW-effect by a factor of ~ 10 for the combination of PLANCK 
with EUCLID, which reaches ~ 7<r. 

(iii) If constraints on cosmological parameters are derived from 
the linear iSW-effect and if contributions from the RS-effect are ne- 
glected, the induced parameter estimation biases are smaller than 
the statistical errors by one order of magnitude because of the 
smallness of the RS-effect in comparison to cosmic variance in- 
duced into the measurement by primary CMB fluctuations and be- 
cause of the strong prior used. If the iSW-effect is used to constrain 
individual cosmological parameters without using a CMB-prior, a 
similar result still applies. Therefore, the RS-effect is negligible as 
a systematic in comparison to other systematics that have been dis- 
cussed in the literature and which have a more pronounced effect 
on cosmogical parameters , e.g. redshif t errors due to peculiar mo- 
tion of the tracer galaxies (Rassa ti2009T) . wea k lensing on the trace r 
population and galaxy magnification bias (Loverde et al. 2007), 
bias evolution of the tracer population dRaccanelli et al.l 1200 8t 
ISchafer et alj|2009h . and contributions due to the kinetic Sunyaev- 
Zel'd ovich effect from reionisation (Gi annantonio & Crittender] 
l2007h . 

Given the estimates that both the spectrum and the bispectrum of 
the RS-effect are only detectable with significances of ~ 0.8cr casts 
doubt on the detectability of this effect in a statistical way, and em- 
phasises the importance of a lternative approaches such as stacking 
methods dGranett et aLll2008h . 
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Figure Al. Nonlinear corrections to the linear CDM spectrum fll \k) (thin 
solid line) of the density field in perturbation theory, at the current cosmic 
epoch: PfP(k) (thin dash-dotted line), the modulus of P[ [ P(k) (thin dotted 



line), and their sum Pss(k) = Pgg\k) 



■ P?P(k) + 2P ( "\k) (thick dashed 



line). The fit to the nonlinearly evolved spectrum Pss(k) from fi-body data 
proposed by Smith et al. 1 2003) is given in comparison (thick solid line). 
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Figure A2. Angular iSW-cross spectrum C Ty (l) of the iSW-effect (solid 
line), split up into the linear effect C^ y \t) (dashed line) and the nonlinear 

(22) ( 1 3) 

RS-corrections C Ty (() + C Ty (I) (dash-dotted line) in logarithmic repre- 
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APPENDIX A: NONLINEAR CORRECTIONS IN 
PERTURBATION THEORY 

Fig. IA1I illustrates the validity of the perturbative corrections to 
P(k) due to nonlinear growt h, by comparison to the result from n- 
body data jSmith et al.ll2003h . Third order perturbation theory is 
able to describe the increase in fluctuation amplitude due to non- 
linear structure formation down to very small scales. One notices 
a deviation between the n-body result and the perturbation theory 
amounting to about 20% in the transition region at a few inverse 
Mpc, corresponding to angular scales of t ^ 300, if most of the 
iSW-signal in the cross-correlation function arises at a comoving 
redshift of x = 2 Gpc/fr, i.e. the maximum of the redshift dis- 
tribution p(z)dz used in this work. The higher orders beyond 3 in 
perturbation theory would correct the difference to the n-body re- 
sult, and it sho uld be kept in mind that the simulation on which the 
description by ISmith et alj d2003l) is based uses slightly different 
cosmological parameters, most notably higher Sl„, and <r 8 . 

The remarkable behaviour of the nonlinear effect to cause an 
anticorrelation between the CMB and the tracer density is shown 
again in Fig. lA2l in a logarithmic representation. 



